Three-dimensional finite element simulation of fracture propagation in rock specimens with pre-existing fissure(s) under compression and their strength analysis

A FORTRAN program, consistent with the commercially available finite element (FE) code ABAQUS, is developed based on a three-dimensional (3D) linear elastic brittle damage constitutive model with two damage criteria. To consider the heterogeneity of rock, the developed FORTRAN program is used to set the stiffness and strength properties of each element of the FE model following a Weibull distribution function. The reliability of the program is assessed against available experimental results for granite cylindrical specimens with a throughgoing, flat and inclined fissure. The calibration procedure of the material parameters is explained in detail, and it is shown that the compressive to tensile strength ratio can have a substantial influence on the failure response of the specimens. Numerical simulations are conducted for models with different levels of heterogeneity. The results show a smaller load bearing capacity for models with less homogeneity, representing gradual coalescence of fully damaged elements forming throughout the models during loading. The maximum load bearing capacity is studied for various combinations of inclination angles of two centrally aligned, throughgoing and flat fissures of equal length embedded in cylindrical models under uniaxial and multiaxial loading conditions. The key role of the compressive to tensile strength ratio is highlighted by repeating certain simulations with a lower compressive to tensile strength ratio. It is proven that the peak loads of the rock models with sufficiently small compressive to tensile strength ratios containing two throughgoing fissures of equal length are similar, provided that the minimum inclination angles of the models are the same. The results are presented and discussed with respect to the existing experimental findings in the literature, suggesting that the numerical model applied in this study can provide useful insight into the failure behaviour of rock-like materials.     

Implementation of a coupled plastic damage distinct lattice spring model for dynamic crack propagation in geomaterials

This paper studies dynamic crack propagation by employing the distinct lattice spring model (DLSM) and 3‐dimensional (3D) printing technique. A damage‐plasticity model was developed and implemented in a 2D DLSM. Applicability of the damage‐plasticity DLSM was verified against analytical elastic solutions and experimental results for crack propagation. As a physical analogy, dynamic fracturing tests were conducted on 3D printed specimens using the split Hopkinson pressure bar. The dynamic stress intensity factors were recorded, and crack paths were captured by a high‐speed camera. A parametric study was conducted to find the influences of the parameters on cracking behaviors, including initial and peak fracture toughness, crack speed, and crack patterns. Finally, selection of parameters for the damage‐plasticity model was determined through the comparison of numerical predictions and the experimentally observed cracking features.     

深水铺管船托管架疲劳分析及动态监测点选取研究

作为S型铺管作业的关键性装备,托管架长期承受着交变载荷的作用,随着海洋开发向超深水发展,结构疲劳破坏问题不容忽视。采用疲劳谱分析的方法并结合线性累积损伤理论,对托管架结构频域下的变形进行了分析,计算了正常海况下和极端海况下托管架疲劳损伤度并对疲劳寿命进行了评估。研究发现托管架在正常海况下作业符合安全要求,在极端海况下局部结构会受到破坏。并从托管架结构安全监测角度,筛选了结构疲劳分析关键点位,为监测点位的选取提供了依据。     

Micromechanical approach to effective viscoelastic properties of micro‐fractured geomaterials

The aim of this paper is to formulate a micromechanics‐based approach to non‐aging viscoelastic behavior of materials with randomly distributed micro‐fractures. Unlike cracks, fractures are discontinuities that are able to transfer stresses and can therefore be regarded from a mechanical viewpoint as interfaces endowed with a specific behavior under normal and shear loading. Making use of the elastic‐viscoelastic correspondence principle together with a Mori‐Tanka homogenization scheme, the effective viscoelastic behavior is assessed from properties of the material constituents and damage parameters related to density and size of fractures. It is notably shown that the homogenized behavior thus formulated can be described in most cases by means of a generalized Maxwell rheological model. For practical implementation in structural analyses, an approximate model for the isotropic homogenized fractured medium is formulated within the class of Burger models. Although the approximation is basically developed for short‐term and long‐term behaviors, numerical applications indicate that the approximate Burger model accurately reproduce the homogenized viscoelastic behavior also in the transient conditions.     

Smeared crack approach: back to the original track

This paper briefly reviews the formulations used over the last 40 years for the solution of problems involving tensile cracking, with both the discrete and the smeared crack approaches. The paper focuses on the smeared approach, identifying as its main drawbacks the observed mesh‐size and mesh‐bias spurious dependence when the method is applied ‘straightly’. A simple isotropic local damage constitutive model is considered, and the (exponential) softening modulus is regularized according to the material fracture energy and the element size. The continuum and discrete mechanical problems corresponding to both the weak discontinuity (smeared cracks) and the strong discontinuity (discrete cracks) approaches are analysed and the question of propagation of the strain localization band (crack) is identified as the main difficulty to be overcome in the numerical procedure. A tracking technique is used to ensure stability of the solution, attaining the necessary convergence properties of the corresponding discrete finite element formulation. Numerical examples show that the formulation derived is stable and remarkably robust. As a consequence, the results obtained do not suffer from spurious mesh‐size or mesh‐bias dependence, comparing very favourably with those obtained with other fracture and continuum mechanics approaches. Copyright © 2006 John Wiley & Sons, Ltd.     

GRADATION AND ASSESSMENT OF DAMAGE DEGREE OF NATURAL LANDSCAPE IN XING, AN COUNTY, GUANGXI

A new technique designed to help quantify the degree of damage to the landscape from one area to another shows a close relationship between population density and the degree of landscape damage. The technique establishes a scale of damage from 0 to 5 (zero = no damage; 5 = severe damage) using data from aerial photographs, land-use maps, and field data. The related formula allows one to compare the relative degree of damage across regions using a combination of an absolute index, a theoretical index, a relative index, and population density. Xing'an County is used to demonstrate the technique.     

Damage identification method using harmonic excitation force considering both modelling and measurement errors

A method of structural damage identification using harmonic excitation force is presented. It considers the effects of both measurement and modelling errors in the baseline finite element model. Damage that accompanies changes in structural parameters can be estimated for a damaged structure from the change between measured vibration responses and ones calculated from the analytical model of the intact structure. In practice, modelling errors exist in the analytical model due to material and geometric uncertainties and a reduction in the degrees of freedom as well as measurement errors, making identification difficult. To surmount these problems, bootstrap hypothesis testing, which enables statistical judgment without information about these errors, was introduced. The method was validated by numerical simulation using a three‐dimensional frame structure and real vibration data for a three‐storey steel frame structure. Copyright © 2005 John Wiley & Sons, Ltd.     

APPLICATION OF MAP AND FILE INFORMATIONVISUALIZATION SYSTEM TO COMPREHENSIVEDIVISION OF NATURAL DISASTERS- Takingthe ChangjiangValleg as an Example

In recent years，scholars at home and abroad have method，and that vanous disasters would be analyzed assynthetically studied natural dlsaste‘theoretically and an Integral．SHI Peilun（1991）putforwad a scientificmethodologlcally，as well as its cases analys。s（CND，term—regional disaster system，which Indicates that1987； MA et al，1990； MARBLE，1990； NE et al，the situation of a disaster（calamity loss）results from1999； PATAK et al，1982； SHI，1991；VAN et al，h…     

Numerical prediction of blast‐induced stress wave from large‐scale underground explosion

This paper presents a numerical model for predicting the dynamic response of rock mass subjected to large‐scale underground explosion. The model is calibrated against data obtained from large‐scale field tests. The Hugoniot equation of state for rock mass is adopted to calculate the pressure as a function of mass density. A piecewise linear Drucker–Prager strength criterion including the strain rate effect is employed to model the rock mass behaviour subjected to blast loading. A double scalar damage model accounting for both the compression and tension damage is introduced to simulate the damage zone around the charge chamber caused by blast loading. The model is incorporated into Autodyn3D through its user subroutines. The numerical model is then used to predict the dynamic response of rock mass, in terms of the peak particle velocity (PPV) and peak particle acceleration (PPA) attenuation laws, the damage zone, the particle velocity time histories and their frequency contents for large‐scale underground explosion tests. The computed results are found in good agreement with the field measured data; hence, the proposed model is proven to be adequate for simulating the dynamic response of rock mass subjected to large‐scale underground explosion. Extended numerical analyses indicate that, apart from the charge loading density, the stress wave intensity is also affected, but to a lesser extent, by the charge weight and the charge chamber geometry for large‐scale underground explosions. Copyright © 2004 John Wiley & Sons, Ltd.     

Simulation of double cold cores of the 35°N section in the Yellow Sea with a wave-tide-circulation coupled model

Based on the MASNUM wave-tide-circulation coupled numerical model, the temperature structure along 35°N in the Yellow Sea was simulated and compared with the observations. One of the notable features of the temperature structure along 35°N section is the double cold cores phenomena during spring and summer. The double cold cores refer to the two cold water centers located near 122°E and 125°E from the depth of 30m to bottom. The formation, maintenance and disappearance of the double cold cores are discussed. At least two reasons make the temperature in the center (near 123°E) of the section higher than that near the west and east shores in winter. One reason is that the water there is deeper than the west and east sides so its heat content is higher. The other is invasion of the warm water brought by the Yellow Sea Warm Current (YSWC) during winter. This temperature pattern of the lower layer (from 30m to bottom) is maintained through spring and summer when the upper layer (0 to 30m) is heated and strong thermocline is formed. Large zonal span of the 35°N section (about 600 km) makes the cold cores have more opportunity to survive. The double cold cores phenomena disappears in early autumn when the west cold core vanishes first with the dropping of the thermocline position. Supported by the National Basic Research Program of China (No. G1999043809) and the National Science Foundation of China (No. 49736190).